Quantum phase transitions of two-species bosons in square lattice 
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We investigate various quantum phase transitions of attractive two-species bosons in a square 
lattice. Using the algorithm based on the tensor product states, the phase boundaries of the pair 
superfiuid states with nonzero pair condensate density and vanishing atomic condensate density are 
determined. Various quantum phase transitions across the phase boundaries are characterized. Our 
work thus provides guides to the experimental search of the pair superfiuid phase in lattice boson 
systems. 



Ultracold atoms in optical lattices provide an unprece- 
dented toolbox for the experimental realization of exotic 
quantum phases which are not easily accessible in con- 
densed matter systems [l|. Recent experimental demon- 
stration of the superfluid-to-Mott transition for ultracold 
bosons in an optical lattice [2[, for example, has paved 
the way for studying other strongly correlated phases 
in various lattice models. Among these systems, Bose- 
Bose, Fermi-Fermi, and Bose-Fermi mixtures have at- 
tracted considerable attention, since a variety of inter- 
esting ordered states can be realized in such systems [l[ . 
In particular, it has been suggested that, for the mix- 
ture of two bosonic species in a lattice, interspecies re- 
pulsion and attraction can give rise to additional in- 
compressible super-counterfluid (SCF) and compressible 
pair superfiuid (PSF) phases, respectively [3-ll7j. Fur- 
thermore, when the intra-species nearest-neighbor repul- 
sion is included, the counterflow supersolid and the pair- 
supersolid phases may also emerge (l7| . 

The quantum phase of our interest is the PSF phase 
of two-species hard-core bosons on a square lattice with 
interspecies attraction U a b(< 0), which is described by 
the Hamiltonian, 
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where a — a,b indicates the two species and n\ is the 
particle densities at site i. The intra-species tunnel- 
ing is denoted by t > and \i is the chemical poten- 
tial. The symbol (ij) indicates nearest neighbors. When 
\U a b\/t 3> 1, the PSF phase can be stabilized by form- 
ing a superfiuid of boson pairs of different species [§-[l3] ■ 
This phase is characterized by a nonzero pair condensate 
density n PSF = \(ab)\ 2 ^ and a vanishing atomic con- 
densate densities Uq = \(cr)\ 2 = 0. By considering two 
species as two components of spin-1/2 particles, it im- 
plies from the symmetry breaking perspective that the 
"charge" U(l) symmetry a — > e l8 a; b — >• e lS b is broken 
spontaneously in the PSF phase, while the "spin" {/(!) 
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symmetry a — »■ e lip a; b — > e~ ttp b is not. Here 6 and ip 
denote the angles of the corresponding {7(1) rotations. 
The PSF phase hence can also be identified by a nonzero 
pair-superfluid density SF oc {(W a + Wb) 2 ) ^ and a 
zero super-counterfluid density pf CF oc {(W a —Wb) 2 ) = 0, 
where W a are the winding numbers for bosons of species 
a [l8|. By contrast, the two miscible superfiuid (2SF) 
phase is expected in the opposite limit of \U a b\/t <C 1, 
where both species form superfluids with nonzero atomic 
condensate density Uq, n h Q ^ 0. In addition to the dis- 
tinct properties of the PSF states, intriguing quantum 
phase transitions out of this phase can occur. Based on 
the analysis of the mean-field theory, it was proposed 
that [l(| , for the cases of spatial dimensions being larger 
than one, the PSF-2SF transitions can be either second 
or first order, while the PSF to Mott insulating (MI) 
phase transitions are always second order. 

While the existence of a PSF phase of lattice boson 
has been pointed out in the literature using various ap- 
proaches [9l-fl7| , direct calculations of microscopic quan- 
tum models with large sizes has not yet been reported ex- 
cept for the one-dimensional case [12l4l4| . This is perhaps 
due to the difficulty in analyzing these quantum models 
by applying conventional methods. For example, despite 
quantum Monte Carlo (QMC) simulation being free from 
the sign problem for this model, several practical techni- 
cal issues may still hinder high-precision large-scale QMC 
calculations. Firstly, multi-worm algorithm needs to be 
implemented 0, El [13; because standard single- worm 
algorithm may converge quite slowly in the PSF states. 
Secondly, very low temperatures are necessary to detect 
the pair superfiuid coherence in the PSF phase. This is 
due to that the effective hopping of pairs is about the or- 
der of t 2 /\U ab \ SI and is quite small since t/\U a b\ <S 1 
in the PSF phase. Finally, systems of large enough sizes 
should be considered in order to determine phase transi- 
tions being of either weakly first order or second order. 
Therefore, accurate quantitative predictions of associated 
quantum phase transitions in the thermodynamical limit 
around the PSF phase remain under investigation. 

In this paper, various quantum phase transitions of the 
model in Eq. ([T]) are investigated. We would set \U a b\ = 1 
as the energy unit. Due to the particle- hole symmetry in 
the present hard-core model, the phase diagram is sym- 
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FIG. 1: (Color online) Ground-state phase diagram of attrac- 
tive two-species hard-core bosons on a square lattice of size 
2 7 x 2 7 as described in Eq. Jl} with \U a b\ = 1- Solid lines 
are guide to eyes. The green dashed line gives the analytic 
estimate of the MI-PSF phase boundary (see text). The inset 
shows the phase diagram in a larger scale. 



metric with respect to the line u = —1/2, at which the 
average density for each species becomes 1/2. It is there- 
fore sufficient to concentrate on the parameter region of 
/j < —1/2 where both species are below half filling. Our 
main purpose is to determine quantitatively the phase 
boundaries of the PSF phase and to reveal the nature 
of these quantum phase transitions. Here the combined 
algorithm [20] of the infinite time-evolving block decima- 
tion (iTEBD) method |2lf and the tensor renormalization 
group (TRG) approach [22| is exploited. The success of 
this method in obtaining accurate physical quantities at 
zero temperature of systems with large sizes has been 
demonstrated for several quantum spin systems [2(| HH- 
[26j . Since systems of two-species hard-core bosons can 
be mapped onto anisotropic spin-1/2 bilayer spin models, 
from our previous experience on the study of the frus- 
trated bilayer spin models (24|, the employed approach 
is expected to give reliable results for the present issues. 
The resulting phase diagram is shown in Fig. [1] The re- 
gion of the PSF phase with finite pair condensate but zero 
atomic condensate is discovered, and the phase bound- 
aries out of the PSF state are determined precisely. The 
nature of various transitions is found to be in agreement 
with that proposed in Ref. [To| . Our findings hence pro- 
vide further theoretical support on searching the PSF 
phases in the bosonic systems with mutual attraction. 
However, we stress that the PSF phase is stabilized only 
within a limited region of the zero-temperature phase 
diagram in Fig. Q] Therefore, carefully tuning system 
parameters into the suggested parameter regime is nec- 
essary to uncover experimentally this novel phase. 

To characterize different phases, several local order pa- 
rameters are evaluated. Because the hopping parame- 
ter and the chemical potential in Eq. ([I} are set to be 
equal for both species, equal densities for the two species 
(n a — n b ) are expected. In our calculations we do find 



that the exchanging symmetry between two species is al- 
ways unbroken and the expectation values for a and b 
bosons are the same. Consequently, we show only the 
results of species a as representatives. When the average 
density of single species boson n — (1/iVs) be- 
comes an integer (here N s denotes the number of lattice 
sites), systems are in the MI states. Otherwise, systems 
belong to either the PSF or the 2SF phases. As men- 
tioned before, the last two phases can be distinguished by 
the expectation values of the atomic and the pair conden- 
sate density n — \(a)\ 2 and tiq SF — |(a&)| 2 , respectively. 
Here (a) = (1/N S ) and (ab) = (1/N S ) In 

the 2SF phase, one has no ^ 0, while in the PSF phase 
one has no = and tiq SF ^ 0. 

These quantities are calculated under the combined 
iTEBD and TRG algorithm [20|, where the ground- 
state wave function is assumed in the form of the ten- 
sor product state (TPS) or the projected entangled-pair 
state (PEPS) [27]. Due to the similarity between the 
present systems and the spin-1/2 bilayer spin models, 
the same construction of TPS as that discussed previ- 
ously in Ref. [24[ for the study of the bilayer systems 
is exploited here. Details can be found in Refs. [23l. [24]]. 
We note that the approximation in TPS can be systemat- 
ically improved simply by increasing the bond dimension 
D of the underlying tensors. The accuracy in evaluating 
the expectation values for a TPS/PEPS ground state of 
very large systems can be systematically improved by in- 
creasing the TRG cutoff Z? C ut- In this study, we consider 
the bond dimension up to D — 5 and keep D cut > D 2 to 
ensure the accuracy of the TRG calculation. As shown in 
the insets of the upper panels of Figs. [2] and [31 the results 
for D = 4 seem to converge to those for larger D = 5. In 
the following, we take D = 4 and D cut = 25 for systems 
with size N s = 2 7 x 2 7 unless mentioned otherwise. 

The results of the local order parameters defined above 
as functions of hopping parameter t for chemical poten- 
tials u = —0.5 and —0.57 are shown in Fig. [2] It is 
found that, at half filling (n — 1/2) with fi = —0.5, 
there exists a second-order quantum phase transition 
at t = < c ,psf-2SF — 0.081 between the PSF and the 
2SF phases. We note that, at half-filling, the particle- 
hole transformation for a single species (say, bi — > b\) 
changes the sign of interspecies interaction in Eq. ([T}, 
U a b — > —U a b- Therefore, at this particular filling, the 
results obtained for the attractive U a b model can be 
reinterpreted in the context of the repulsive U a b model 
and vice versa. In particular, the PSF order parameter 
(ab) is mapped onto the super-counterflow (SCF) order 
parameter (aw). Therefore, the critical value of t for 
the PSF-2SF transition of the attractive hard-core bo- 
son model should be the same as the SCF-2SF transi- 
tion of the repulsive hard-core boson model studied pre- 
viously [H 0] • We find that our critical value of t is con- 
sistent with the result (i c ,psF-2SF — 0.092) obtained by 
QMC method Q. This indicates that the bond dimen- 
sion used here is large enough to provide accurate find- 
ings for the present model. As shown in Fig.^a), within 
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FIG. 2: (Color online) Values of particle density n, atomic 
condensate density no, and pair condensate density tiq SF for 
the ground states at (a) u = —0.5 and (b) —0.57 as functions 
of hopping parameter t. The inset in (a) illustrates the D 
dependence of no, where we use D cut = 16 for D = 3 and 
Dcut = 25 for D = 4, 5. 
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FIG. 3: (Color online) Upper panel: values of n, no and no SF 
for the ground states at t — 0.12 as functions of chemical 
potential fi. The inset shows the D dependence of no, where 
we use Dcut = 16 for D = 3 and D cut — 25 for D = 4, 
5. Lower panel: n and no as function of u at different t's. 
The arrows denote the analytic predictions of the transition 
points, assuming that they are continuous transitions. 



the PSF phase, the values of the pair condensate density 
tIq SF do not remain the same and are somewhat smaller 
than the constant value n FSF = 1/4 predicted by the 
mean- field theory This reduction should result from 
quantum fluctuations therein. 

While the PSF-2SF transition at fj, = —0.5 is of second 
order, based on the mean-field analysis, it was suggested 
that this transition could become first-order around the 
end point of the PSF-2SF phase boundary If this 

conclusion is true, the behavior of the first-order tran- 
sition should become more apparent as the end point 
is approached. In Fig. Hfb), n, no, and n FSF as func- 
tion of t for fi = —0.57 are displayed. A transition at 
£ Cj psf-2SF — 0.1024 is found, beyond which no becomes 
nonzero. Besides, anomalies in n and n FSF occur at this 
transition point. Although this transition point is quite 
close to the end point of the PSF-2SF phase boundary 
as indicated in Fig. [TJ we find that the transition is of 
second order. The same conclusions about the order of 
transition are reached for larger /x's. Nevertheless, we 
do observe that the order-parameter curves become more 
and more first-order-like as /j, decreases from /i = —0.5 to 
u = —0.57 (not shown here). Our observation is consis- 
tent with the results achieved by the Monte Carlo simula- 
tion on the related classical model in Ref. [l(| , where only 
the second-order PSF-2SF transitions are found. How- 
ever, one cannot completely rule out the possibility of 
weakly first-order transitions with tiny jumps in order 



parameters when the transition points get much closer 
to the end point of the PSF-2SF phase boundary. 

As seen in Fig. [D^b), when fi = —0.57, an additional 
second-order phase transition between the n = MI 
and the PSF states appears at i c ,Mi-PSF — 0.094. In 
this case, the insulating state is expected to be unsta- 
ble when the energy gap of pair excitation closes as t 
increases. Above the transition point, the lowest pair 
excitation starts to condense and a PSF state devel- 
ops. By using the expression of the two-particle exci- 
tation energies relative to the n = state up to the 
third-order perturbation expansion in t [l6l |28| , the crit- 
ical hopping parameter for this transition is found to be 
ie,Mi-PSF = v ( — 2/x — l)/4z, where z — 4 is the coordi- 
nation number for square lattices. As shown in Fig. HJTd) 
and Fig.[T]for cases of other /i's, our findings of £ c ,mi-psf 
agree well with the values given by the analytical ex- 
pression. This further substantiate that the employed 
method can provide reliable results for the present sys- 
tem. 

When the single-particle hopping parameter t is suffi- 
ciently large, instead of the MI-PSF transition, a direct 
transition from the n = MI state to the 2SF state 
as chemical potential \x being increased becomes pos- 
sible. If this MI-2SF transition is of second order, it 
should be induced by closing the energy gap of single- 
particle excitation. By using the analytical expression of 
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the one-particle excitation energies relative to the n — 
state [lg, [28[ , the critical value of n at fixed t is given 
by A* c ,mi-2SF = —zt. However, as discussed in Ref. fioj . 
the continuous MI-2SF transitions can be preempted by 
the first-order ones when mutual interaction of the di- 
lute gapped quasiparticle excitations in the MI phase is 
attractive. Due to the energy gain from the mutual at- 
traction among quasiparticles, this first-order transition 
can happen before closing one-particle excitation gap. 
That is, the transition points for the first-order transi- 
tions would be lower than those derived by the gap clos- 
ing condition. These discussions are confirmed by our 
calculations. In Fig. [3l our results of the local order pa- 
rameters n, no and rig as functions of a for various t 
are plotted. For t = 0.12 (upper panel), a first-order MI- 
2SF transition at /x c ,mi-2SF — —0.611 is clearly observed. 
Upon increasing t, the jumps in the order parameters be- 
come smaller and smaller, and the transition eventually 
turns to be of second order (see the lower panel). For 
example, the MI-2SF transition at t = 0.6 is found to be 
continuous, and it occurs at /i c ,Mi-2SF — —2.4 in agree- 
ment with the analytic formula. As shown in Fig. [31 
the transition points for those first-order transitions are 



always below the values given by the gap closing condi- 
tion. Thus the phenomena of the transitions preempted 
by first-order ones is indeed demonstrated in the present 
systems. 

To conclude, various quantum phase transitions of the 
two-species hard-core boson model with attractive in- 
terspecies interaction are explored under the combined 
iTEBD and TRG algorithm (20|. Clear evidence of a 
PSF phase over a finite regime in the phase diagram is 
provided. Critical values of system parameters at cor- 
responding first- or second-order phase transitions are 
evaluated. Besides providing quantitative predictions of 
the phase boundaries of the proposed PSF phase, the 
present work also illustrates clearly the general validity 
and flexibility in applying the current formalism to de- 
termine quantum phase transitions in two dimensions. 
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